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1. Introduction 

Over the last 30 years there has been an increasing interest in the critical behaviour of 
systems (fluids, magnets . . . ) confined between two infinite parallel plates i.e. films. Such 
systems can be regarded as (i-dimensional generalization of two {d — l)-dimensional walls 
separated by matter of thickness L. The so called thermodynamic Casimir effect related to 
fluctuation induced long-ranged forces between the plates was, and still is, the central theme 
of these investigations both experimentally and theoreticall}j§|. This effect was first predicted 
by Fisher and de Gennes (Fisher & de Gennes 1978) in their investigation on a confined 
critical binary liquid at its bulk demixing point. They showed that the reduced free energy 
per unit area contains a term of the form A'^^^L"'^"'"^, where A*^^) is the so called Casimir 
amplitude by analogy with the Casimir effect in vacuum fluctuations of the electromagnetic 
field between two metallic plates (Casimir 1948). It is a universal quantity that depends on 
the bulk universality class and the boundary conditions (r) imposed on the confining walls 
(Privman & Fisher 1984, Singh & Pathria 1985a). The different types of boundary conditions 
are in turn related to distinct universality classes of surface critical behaviour depending on 
the behaviour of the order parameter at the surfaces bounding the system and some additional 
surface properties (Diehl 1997). 

In general, according to Privman & Fisher (1984) and Singh & Pathria (1985a) the 
singular part of the free-energy density of a finite (in one or more directions) rf-dimensional 
system with linear size L, near the bulk critical point T^, may be expressed in the form 

fi;}{t,h-L) ^ L-'Y^^\c,tL'''',C2hL^'^), (1.1) 

where t and h are related to the temperature, T, and the external magnetic field, H, via 

t = h = j^. (1.2) 

The arguments of Y^'^\xi,X2) are appropriate scaled variables, v and A are the usual critical 
exponents, while all the details of the system are incorporated in the non-universal quantities 
ci and C2- Then the function Y^'^\xi,X2) is a universal scaling function, whose exact 
expression depends upon the number of finite directions, the bulk universality class and the 
boundary conditions (r) to which the system is subjected. For a system confined to a film 
geometry, at T = T^, the Casimir amplitude coincides with the critical amplitude of the 
singular part of the free energy density, i.e. 

A(") = F(")(0,0). (1.3) 

The spherical model of Berlin and Kac (Berlin & Kac 1952) was initially designed to 
mimic the critical properties of the Ising model. It has been obtained by requiring the spins 
to be continuous variables subject to a global constraint (the sum of the squares of spins at 
each lattice site is equal to the total number of sites N) rather than a local one (the square of 
the spin at each site is exactly 1). Later it was shown that the free energy of this model can 

§ The literature on the thermodynamic Casimir effect can be found in (Griineberg & Diehl 2008) and references 
therein. 
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be obtained as a limiting case of that of the Heisenberg model with infinite number of spin 
components (Stanley 1968). The equivalence between the Heisenberg model with infinite spin 
components and the spherical model remains valid for finite systems, as long as one considers 
boundary conditions that preserve the translational invariance of the lattice (Knops 1973). 
Lewis & Wannier (1952) proposed to simplify the spherical model by requiring the global 
constraint of Berlin and Kac to be obeyed in the sense of an ensemble average. This model is 
generally known as the mean spherical model. 

Because of its exact solubility the ferromagnetic mean spherical model has been 
extensively used to gain insights into the critical properties of finite (in one or more directions) 
system^. Periodic boundary conditions, implying ferromagnetic interactions of spins 
belonging to both boundaries on the finite directions, have been by far the most used boundary 
conditions. These allow for analytic treatment of problems related to the finite-size scaling 
theory (Privman et al. 1991, Brankov et al. 2000). In addition to periodic boundary conditions, 
antiperiodic boundary conditions have also been used (Barber & Fisher 1973, Singh & 
Pathria 1985^?) to investigate the finite-size scaling properties of the ferromagnetic mean 
spherical model. The aforementioned boundary conditions do not break the translational 
invariance of the model in the absence of a magnetic field. Much less has been achieved in 
the case of nonperiodic (free) boundary conditions. These are believed to be more relevant to 
real systems, especially to systems confined between parallel plates. The ferromagnetic mean 
spherical model of finite thickness has been investigated in the case of Dirichlet boundary 
conditions by Barber & Fisher (1973), Barber (1974), Danchev et al. (1997), Chen & Dohm 
(2003) and Dantchev & Brankov (2003), and Neumann boundary conditions by Barber et al. 
(1974), Danchev et al. (1997) and Dantchev & Brankov (2003). Barber & Fisher (1973) 
and Barber (1974) investigated the scaling properties of the mean spherical model with a 
finite thickness using a method originally devised by Barber & Fisher (1973). They obtained 
explicit forms of the equation for the spherical field at three, four and five dimensions, 
separately. The method was extended to the study of the finite-size effects in the case of 
Neumann and Neumann-Dirichlet boundary conditions by Danchev et al. (1997). Chen & 
Dohm (2003) argued that the results of Barber & Fisher (1973) and Barber (1974) for Dirichlet 
boundary conditions were incorrect at three dimensions far from the critical point. Later, 
the equation for the spherical field of Barber & Fisher (1973) was rederived by Dantchev & 
Brankov (2003). It should be mentioned, however that the derivation of Dantchev & Brankov 
(2003) is based on an improved method of Barber & Fisher (1973). On the other hand, in 
(Barber & Fisher 1973, Barber 1974, Danchev et al. 1997, Dantchev & Brankov 2003) the 
free energy density for the different kinds of boundary conditions was obtained by integrating 
the equation for the spherical field leading to complicated integral representations. 

In the present paper we propose a different method to treat the finite size effects in the 
mean spherical model of finite thickness. The method is a generalization of that devised 
by Singh & Pathria (1985a) for periodic boundary conditions. It applies to antiperiodic, 
Dirichlet, Neumann boundary conditions imposed on both surfaces bounding the system and 

II For an extensive list of literature see (Barber & Fisher 1973, Singh & Pathria 1985fl, Privman 1990, Chamati 
et al. 1998, Brankov et al. 2000) and references therein. 
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a combination of Dirichlet and Neumann boundary conditions on each surface. The method 
is quite general and is used for arbitrary dimension. Furthermore, the finite-size contributions 
to the free energy density are obtained directly form the corresponding general expression 
rather than through integration of the equation for the spherical field. We anticipate here that 
from our results we recover the particular cases of Barber & Fisher (1973), Barber (1974), 
Danchev et al. (1997) and Dantchev & Brankov (2003) but not the results of Chen & Dohm 
(2003). We will return to these points later in the paper. 

The rest of the paper is structured as follows: In Section [2] we define the model and 
present the expressions for the free energy density and the equation for the spherical field 
for the different kinds of boundary conditions. In Section [3] we present in detail how the 
method applies to the case of periodic boundary conditions and compare our results with 
those available in the literature. In Sections |4] through |7] we extend the method for arbitrary 
dimensions to the other boundary conditions and compare with the results obtained by other 
authors using different methods. Finally in Section[8]we summarize our results. 

2. The model 

We consider the mean spherical model on a (i-dimensional lattice confined to a film geometry 
i.e. infinite m d — \ dimensions and of finite thickness L in the remaining dimension with 
volume V = L oo'^^^. The linear size of the lattice is measured in units of the lattice 
constant, which will be taken equal to 1 . The model is defined through 



where Si = Si{r, z), the spin at site i, is a continuous real variable (— oo < Sj < oo) with 
coordinates r in the d — 1 infinite "parallel" planes and z in the finite lateral direction. Jij is 
the interaction matrix between spins at sites i and j, and H is an ordering external magnetic 
field. Finally the field /i is introduced so as to ensure the spherical constraint 



where (• • ■) denotes standard thermodynamic averages computed with the Hamiltonian 7i and 
N the total number of spins on the lattice. 

Along the finite z direction we impose different kinds of boundary conditions which we 
will denote collectively by r. For a lattice system this means: 

[p] periodic: s(r, 1) = s(r, L + 1); 

(a) antiperiperiodic: s(r, 1) = — s(r, L + 1); 

{D) Dirichlet: s(r, 0) = s(r, L + 1) = 0; 

(A^) Neumann: s(r, 0) = s{r, 1) and s{r, L) = s{r, L + 1); 

(ND) Neumann-Dirichlet: s{r, 0) = s{r, 1) and s{r, L + 1) = 0. 




(2.1) 




(2.2) 
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For nearest neighbour interactions and under the above boundary conditions the 
Hamiltonian (12.11) may be diagonalized by plane waves parallel to the confining plates and 
appropriate eigenfunctions ipn\z) using the representation 



dqi f^^ dqa- 



2n 



271 



(2.3) 



where the integrals over g^s, the components of q, are restricted to the first Brillouin zone of 
the hypercubic lattice of dimension d — 1. The orthogonal eigenfunctions read 



1 
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exp 
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(2.4a) 
(2.4Z7) 
(2.4c) 

(2.4J) 

(2.4e) 

(2.5a) 

(2.5Z7) 

(2.5c) 
(2.5J) 
(2.5c) 



Those corresponding to the diagonalized interaction matrix are An — 2 + cjg \ 

At zero field, the free energy density of the mean spherical model under the above 
boundary conditions imposed along the finite lateral size has the general expression (obtained 
via the Legendre transformation) 

(3F^J\t, L;<P) = -\nK + —J2 U,^M + ' 4^^) - ^^i^ + Jj^) - K, (2.6) 



n=0 
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where K = pj = J/keT and 




d 



In z + 2^(1 -cos ^i) 



(2.7) 



i=l 



In (|2.6I) . the shifted spherical field, defined by 



(j) = ji/ J — \ 



is the solution of the spherical constraint 



- LU,, 




(2.8) 



n=0 



with the bulk function 




dOd 1 

2vr^ + 2Eti(l-cos^^.)' 



(2.9) 



whose asymptotic behaviour has been studied in considerable detail for z G C in (Barber & 
Fisher 1973). 

The bulk limit is obtained from (|2.8I) by letting the lateral size L go to infinity. This 
allows us to investigate the thermodynamics of the ferromagnetic mean spherical model in 
the thermodynamic limit. Here we shall not enter into the investigation of the thermodynamic 
properties (the interested reader may refer to (Pathria 1996)). It is worth mentioning that the 
present model undergoes a continuous phase transition at a bulk critical point determined by 



for c? > 2. On the other hand for d > A the model exhibits mean-field-like critical 
behaviour. Recently an efficient method to estimate the Watson integral Wd(0) and the 
associated logarithmic integral for (i-dimensional hypercubic lattice has been proposed by 
Joyce (feZucker (2001). 

Before embarking into the investigation of the finite-size effects for the different 
boundary conditions a few comments are in order: 

(1) For the evaluation of the finite-size contributions to the bulk expressions of the free 
energy and the equation for the spherical field the case of periodic boundary conditions is the 
simplest. A powerful method to treat these sums was proposed in (Singh & Pathria 1985a). It 
was found that in this case the correlation length of the finite system is given by = (f)~^^'^ 
(Singh & Pathria 1987). 

(2) The same method was successfully extended to the case of antiperiodic boundary 
conditions (Singh & Pathria 1985Z7). However, the calculation were done in the absence of 
the lowest mode uj^^^ in (12.61) and (12.81) . Later by investigating the correlation function (Allen 
& Pathria 1993), it was found that the correlation length and the root square of spherical 
field (f) are no more connected by a simple relation as it is the case for periodic boundary 
conditions. It was suggested that the correlation length is in fact related to the solution of the 
equation for spherical field shifted by the asymptotic behaviour of lUq i.e. j^. 



K,,d = W^d(O) 



(2.10) 
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(3) The remaining boundary conditions have attracted less attention (Barber & Fisher 
1973, Barber 1974, Barber et al. 1974, Danchev et al. 1997, Brankov et al. 2000, Chen & 
Dohm 2003, Dantchev & Brankov 2003). Apart from the paper by (Chen & Dohm 2003), all 
investigations were specialized to (i = 3 using an approach based on (Barber & Fisher 1973). 
For arbitrary d and Dirichlet boundary conditions, a recent work by Chen & Dohm (2003) 
proposed a different method and claimed that the results of (Barber & Fisher 1973), obtained 
at d = 3, were incorrect. 

Here we generalize the method of (Singh & Pathria 1985a) for periodic boundary 
conditions to the other boundary conditions. We investigate the finite-size effects of the free 
energy density (12.61) and the equation for the spherical field (12.81) in each case for arbitrary 
dimension and comment on the results of the aforementioned papers. 



3. Periodic boundary conditions 



For the sake of completeness we will derive here the relevant expressions for the 
mean spherical model of finite thickness under periodic boundary conditions for arbitrary 
dimensionality. The derivation is adapted from (Singh & Pathria 1985a) (see also (Chamati 
et al. 1998)). The explicit expressions for the free energy density (|2.6I) and the equation for 
the spherical field (|2.8I) are given by 

'27r 



n=0 



1 — COS 



L 



-n 



and 



K 



1 



n=0 



1 — COS 



2tt 



-n 



(3.1a) 



OAb) 



respectively. In this case uo^^ = 0. 
Using the integral representations 



-zt 



and 



In z 



/"°° df 



(3.2a) 



(3.2b) 



we can separate the expressions for the free energy density (I3.1al) and the equation for the 
spherical field (I3.1Z?I) into the corresponding bulk expressions and the associated finite-size 
contributions. Let us illustrate how this works for the equation for the spherical field. With 
the aid of (|3.2al) the sum entering this equation can be written as: 

L-l 



s. 



n=0 

2" dOi 

27 



2n 

+ 2 I 1 — cos —n 



2n 



d9^ 
2tx 



dz exp < —z 



+ 2 + 2^(1 -cos^, 
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with 



L-l 



ra=0 



27r 
z COS —n 



(3.4) 



Notice that the summand in (13.41 ) is a periodic function of period 27r. Further we use (a 
generalization of) the Poisson summation formula, namely 

b oo 

(3.5) 



to get the identity 

Qt')(^) = X;exp 



n=0 



27r 
2; cos —n 
Ij 



(3.6) 



l=—oo 



where luix) stands for the modified Bessel function of the first kind (Abramowitz & 
Stegun 1972). 

Substituting (|3.6I) into (|3.3I) we obtain 

Stl{^)=LYl / dze-^^[e-^-^I,{2z)Y-\''^Iu{2z) 
1- ^ Jo 



l=—oo ' 

= W,(0) +2L V / rf^e-^^ [e-'^Io{2z)Y-\-'^Iu{2z), 
1=1 -^0 

where we have used the integral representation 

POO 

WM= / dze-''^[e-'''Ioi2z)Y 
Jo 

and the relation I-2p{x) = hpi^) (Abramowitz & Stegun 1972). 

For large L, using the asymptotic expansion (Singh & Pathria 1985a) 

x-u^/2x r 1 n QO,,2 



(3.7) 



(3.8) 



1 9 - 32z/^ 
^^8^^ 2!r8x)2 



v2vnr 

after some straightforward steps, keeping only leading terms in we get 

4 



5^5(0) = Liy,(0) + L- 



(3.9) 



(3.10) 



(4vr)t ' tr (l^^v^)^-^ 
where K^{x) is the modified Bessel function of the second kind (Abramowitz & Stegun 1972). 
Collecting the above results the equation for the spherical field reads 

4 



K = Wd{(p) + 



(3.11) 



(4vr)t' tr {¥L^f^-' 
The finite-size behaviour of this equation has been studied in great detail by Chamati et al. 
(1998) for arbitrary d for different geometries including questions like dimensional crossover, 
finite-size shift of the bulk critical temperature, etc. In particular it was found that the finite- 
size shift obeys the predictions of the finite-size scaling. For 2 < d < A, using the asymptotic 
behaviour 



WAz) 



(4vr)d/2- 



2-d 



(3.12) 
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equation (13.1 II) takes the scaling form 



y 



.d-2 



(47r)'^/2 



2-d 



(3.13) 



1=1 \2^yY 

where we have introduced the scaling variable y = L-/(j) and x = L^/^^Kc^^ — K) with 
u = (d — 2)^^ - the critical exponent measuring the divergence of the correlation length. 
Consequently we have a solution of the form = 0^^^^ = Lfp{>c), where /(p) is a universal 
scaling function. For arbitrary d the nature of the scaling function f(p){x) can be determined 
only numerically (see e.g. (Chamati et al. 1998)). For the particular case = 3, equation 
(13.131) takes a simple form 

y 



Itxk = ln2 sinh — . 

2 



(3.14) 

The solution of this equation leads to the universal scaling function 

y = g(^p){x) = 2 arcsinh (le^^^) (3.15) 

At the critical point, x = 0, we obtain the critical amplitude of the finite-size correlation 
length ^L- 

1 + 



yo = 9ip){0) = 21n- 



(3.16) 



The finite-size contributions to the free energy for any d can be accounted for by using 
the integral representation (13. 2M . The aim here is to transform the sum 

L-l 



n=0 
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1^ 27 



n=0 
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(3.17) 



into a more tractable form suitable for analytic treatment. After some straightforward algebra 
along the lines explained above, including the use of the identity (|3.6I) . we arrive at 

V^^l{<P) = LUM)-L^jiY.Yi 



(47r)f 



(3.18) 



1=1 {\IL^Y 

Using (|3.18l) for the large L asymptotic behaviour of the free energy density (I3.1al) we 



get 



(3F^f\T,L- 



where 



F,{T- 0) = lim F^;\t, L;(j)) = — [InK + U,{<P) - K<P - 2K] 



(An 



1=1 {\IL^Y 



(3.19) 



(3.20) 



is the bulk free energy density. 

For 2 < d < A,'m the vicinity of the bulk critical point, we have the expansion 



(47r)^/2 d 



2 - d 



yd/2 



0{z 



(d+l)/2N 



(3.21) 
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Then the singular part of the free energy density (13.191 ) takes the scaling form 



(47r)'^/2 1 d 



2-d 



1=1 



K, {ly) 



(3.22) 



where y is the solution of (13.131) . Thus, the scaling behaviour of the free energy density is 
consistent with the finite-size scaling hypothesis (|l.ll) . The Casimir amplitude for d = 3 
i.e. the critical amplitude of the singular part of the free energy density can be computed 
analytically, having in mind the solution (13.161) . with the aid of polylogarithmic identities 
(Sachdev 1993). The result is (Sachdev 1993, Danchev 1998): 

2C(3) 



A(P) 



Svr 



-0.153051... 



(3.23) 



This value is compatible with estimations obtained for more realistic 0{n) models using 
different approaches. For more details we refer the reader to (Griineberg & Diehl 2008). 

As one would expect from previous studies on the mean spherical model, the results 
obtained for periodic boundary conditions are in conformity with the finite-size scaling 
hypothesis. In the following we will extend the method described here to the other 
boundary conditions. The generalization is somehow straightforward, however we will see 
the appearance of some subtleties that need careful consideration. 



4. Antiperiodic boundary conditions 



The explicit expressions for the free energy density (12.61) and the equation for the spherical 
field (12.81) for arbitrary dimension read 



(T,L; 



1 

2L 



L-1 



d-l 



n=0 



+ 00, 



1 



K 



and 



L-1 



K 



n=0 



respectively. Here the lowest eigenmode cug"'' 7^ 0, in contrast to the case of periodic boundary 
conditions. 



-K, (4.1fl) 



{A.lb) 



Unlike the analysis of (Singh & Pathria 1985Z?) we shall not omit uJ^ 



from our equations, rather we will use the combination a 



(r) 


,(«) 



-2 + 2 cos f 
= + uj^^ as a variable in 
our consideration of the sums entering (|4.1al) and (|4.1Z?I) . We will see that this is crucial to 
our further treatment. Indeed, in our notations, is expected to define the correlation length 
and no other definition for this quantity is necessary as it has been suggested by Allen & 
Pathria (1993), where the behaviour of order parameter correlation function for the finite 
mean-spherical model with antiperiodic boundary conditions has been investigated in detail. 
The analysis of Section [3] for periodic boundary conditions can be applied to the sums 

=E^'^-i[^^'^^-^""^] (4.2a) 

n=0 



Finite-size effects in the spherical model of finite thickness 



11 



L-l 



{A.lb) 



n=0 



appearing in the left hand side of (I4.1al) and (|4.1Zp|) . respectively. Now instead of the identity 
(Il6l) we use (Singh & Pathria \9%5b) 
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n=0 



Z cos ^(■'^ + 2/ 
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L ^ cos{7iI)Ili{z) 



(4.3) 
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to end up with the final expressions 
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(4.4a) 
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(4.4Z7) 



+ ^0 



For 2 < d < 4, making use of the asymptotic expansion (13.121) . for large L and in the 

he scaling behavioi 



d-2 
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(47r)'^/2 





2 - 



vicinity of the bulk critical point, with (j) + lJq < 1, we have the scaling behaviour 



1=1 



-1 



(4.5) 



where we have used oo. 



(a) 



Following the analysis of (Allen & Pathria 1993) it is easy 
to show that coincides with the inverse of the finite-size correlation length ^l. In this 
case, the solution of (14. 5|) may be written as = Lf(^a){^), where /(„) is a universal scaling 
function. 

The critical temperature of the film corresponds to y = n i.e. = (^) , which is the 
asymptotic of the lowest mode cUq""* for large L. Setting y = tt in equation (|4.5I) . we find that 
the critical point of the film is shifted from the bulk one by a quantity proportional to L~^^'^ 
in agreement with the finite-size scaling predictions. 

For arbitrary dimension d, equation (|4.5I) can be solved only numerically. Here we will 
specialize to the three-dimensional system, which allows analytic treatment. For d = 3, 
equation (|4.5I) transforms into 



2tih = In 2 cosh - ^/y^ 



whose positive solution reads 

y = 9{a){^), 



(4.6) 
(4.7) 
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with the universal scaling function 

^?(,)(x) = (vr^ + 4[arccosh {\e^^'')ff'' . (4.8) 
At the bulk critical point, i.e. x = 0, we obtain 

yo = giam = ^n. (4.9) 

This is the critical amplitude of the finite-size correlation length ^l. This result has been 
obtained also in (Allen & Pathria 1993) using a different definition for the correlation length 
imposed by the choice of a different initial equation for the spherical field. 

Using the asymptotic behaviour (13.211) . valid for 2 < (i < 4, and the fact that y is the 
solution of equation (14.51 ). we may write the singular part of the free energy density (|4.1a|) in 
a scaling form as 



1 



+ 



TT 



TC 



(47r)<i/2 



d 



-r 



V 



2-d' 



\l 2 



iVy' 



1=1 



(4.10) 



In accordance with the finite-size scaling hypothesis (11.11) . For the important case d = 3, its 
critical amplitude at the bulk critical temperature and consequently the corresponding Casimir 
amplitude for antiperiodic boundary conditions i^ 

A^"") = 0.274543 ... . (4.11) 

Notice that the Casimir amplitude here is positive in contrast to the case of periodic boundary 
conditions, see e.g. (13.231) . but approximately twice higher in magnitude. This result is 
compatible with that of (Krech & Dietrich 1992) obtained using renormalization group. 

Before closing this section let us mention that the expressions for periodic boundary 
conditions and those corresponding to antiperiodic boundary conditions may be written in 
unified general forms with a parameter characteristic of twisted boundary conditions. In the 
remainder of the paper we will use these methods to analyse the finite-size effects of the mean 
spherical model with a film geometry subject to more realistic boundary conditions. 



5. Dirichlet boundary conditions 

5.1. Equation for the spherical field 



Here we will evaluate the finite-size contributions of the sum appearing in (12.81) in the case of 
Dirichlet boundary conditions for arbitrary d. We start with 

L 



sf^^) = Y.^<i-i 



n=l 



+ 2 COS 



TT 



L + l 



2 cos ■ 



vm 



L + l 



(5.1) 



^ This result was obtained independently by (Dantchev & Griineberg 2008) 
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which upon extending the sum to n = 2L + 1 may be written 

2L+1 



=0 



n=0 

2' 



+ 2 cos 
+ 2 cos 



TT 



— 2 cos 



L + l 



1 

+ 2 



vrn 



+ 2 COS 



TT 



L + l 



(5.2) 



The last two terms correspond to n = L + 1 and n = 0, respectively. Comparing with (13.31) 
we find that the sum in the right hand side is exactly 5^^2^_,_2 (^(p + f^o^^ j corresponding to 
periodic boundary conditions with a film of thickness 2L + 2. This suggests that the analysis 
of the sum appearing in (|5.2I) can be performed following the method outlined in Section [3] 
(see e.g. (ITTOl) ). 

The equation for the spherical field follows from (|2.8I) . (15.21) and (13.101) . This is (to the 
leading order in L^^) 



K = W, 



+ 



TT 
Z2 



1 

L 



W, 



71 
l2 



(4vr)rf/2 



-1 oo 



-1 



2 
2/L 



-W, 



d-l 



l2 



+ 4 



-W, 



d~l 



TT 
l2 



L2 



1=1 



(5.3) 



L2 



where we have used the large asymptotic behaviour uj^^^ ~ — fi-. Equation (15.31) is the general 
form of the equation for the spherical field for arbitrary dimension d. Notice that the right hand 
side is composed of a bulk term, a size dependent surface term 



W, 



1? 



71 
Z2 



4 1 - ^Wd-i 



(5.4) 



and finite-size corrections. Here, all the quantities are function of the combination — 
Thus, the critical properties of the mean spherical model of finite thickness under Dirichlet 

2 

boundary conditions should be investigated using — ^ as a small parameter keeping in mind 
that < 1 and L > 1. 

The scaling behaviour of (15.31) depends strongly upon d, indeed the asymptotic expansion 
of Wd~i{z) for small argument takes different expressions for different values of d. We 
first start with the important three dimensional case that has been extensively studied in the 
literature. Later we will extend our analysis to the interval 3 < < 4. This constraint for d 
ensures the validity of the asymptotic expansion (13.121) for Wd{z) and Wd-i{z) at the same 
time. 

For d = 3, using (13.121) with d = 3 and 

^ (51n2-ln2) +0{z\nz), 



W^iz 
form (15.31) we have 



47r 



1 

L 



47rL 



InL - 



47rL 



In 



sinh2^L20-7r2 



(5.5) 



(5.6) 
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which, apart from an unimportant numerical factor that enters in the definition of K, coincides 
with equation (4.69) of (Barber & Fisher 1973) obtained at d = 3 via a different method. 
Consequently our result disagrees with that of (Chen & Dohm 2003), where it has been argued 
that equation (4.69) of (Barber & Fisher 1973) was incorrect far from -K^cJll. At the bulk 
critical point, solving (|5.6I) . with the assumption La/0 <^ 1, we find (p ~ in agreement 
with (Chen & Dohm 2003, Dantchev & Brankov 2003). It is worth mentioning that (15.61) 
cannot be put in a scaling form because of the logarithmic dependence on L. However, if 
one considers a shifted critical point that absorbs the term proportional to In L and the surface 
contributions according to (Barber & Fisher 1973) 



K 



(D) 
s,3 



Kc,3 + -r 



K, 



c,3 



71n2 



57r 



InL, 



(5.7) 



it is possible to recover the scaling behaviour. Notice that expression (15.71 ) does not conform 
with the predictions of the theory of finite-size scaling for the finite-size shift of the critical 
temperature temperature, since there is a term containing In L. 

In the case 3 < d < A, using the asymptotic expansion (13.121) . the equation for the 
spherical field (15.31) may be written as 

[WM - |W^d-i(0) - |W^d-i(4)] 



+ 



+ 



(4vr)d/2 
1 

(47r)^/2 



r 



2-d 



3-d 



-1 



°° Kd 



71' 



-1 



(5.8) 



where we have introduced the usual notations x = L^/^iyKc^^ — K) and y = Ly/cj). The first 
term in the right hand side shows that the standard finite-size scaling hypothesis breaks down 
in the vicinity of the bulk critical temperature. It is likely that the introduction of a scaling 
function depending on x and an additional variable to take care of the term proportional to 
L'^^^ aiming at the modification of the finite-size scaling would cure this deficiency. 

From equation (15.31) it is easy to see that leading large L behaviour of the finite-size shift 
from the bulk critical temperature to the one where the film is expected to have a singular 
behaviour is L^^. This result, valid for any dimension, shows that the finite-size scaling 
hypothesis is violated as it has been pointed out by Brezin (1983). 



5.2. Free energy density 



Let us now turn to the evaluation of the free energy density (|2.6I) . We need the asymptotic 
behaviour of the sum 

L 



n=l 



+ 2 cos 



TT 



L + l 



2 cos 



7m 



L + l 



(5.9) 



The origin of this discrepancy is discussed in (Chen & Dohm 2003). 
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entering the expression (12.61) of the free energy. Here again we extend the sum to n = 2L + 1, 
to get 

H 2L+1 r 

vr _ Tin 



n=0 



+ 2 COS 



L + l 



— 2 cos 



L + l 



+ 2 cos 



TT 



L + l 



2 cos ■ 



TT 



- 2 



,(p) 



The sum in the right hand side of (15.101) corresponds to 2L+2 
Then from (13.181) and (15.101) we get the free energy density 



4- i 



(5.10) 
from (IOtI) . 
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djSurf . 



I2 



(4vr) 



Z2 



isTd 2/L 



(5.11) 



where 



1 



1 



f/d(^)-7^f/d-l(^ + 4)--f/rf_i((T) 



2 - -V / 2 

accounts for size dependent contributions stemming from the surfaces. In the following 
we will discuss the finite-size behaviour of the free energy and its dependence on the 
dimensionality d. Again we discuss separately the cases d = 3 and 3 < d < 4 imposed 
by the validity of the expansion (13.211) for \Jd[_z) and \J^^\[^z), simultaneously. 
For d = 3, we use the expansions (|3.21l) and 



\J2{z) = f/2(0) - ■^z\nz + ^{l + 5ln2)z + 0{zHnz) 

4:71 4:71 



to obtain the singular part of the free energy density 

3 1 



-'- / 2 2\ 1 

-X (y — vr 

2 ^ 12vr 
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16vr 
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' 8vr 



vr 



- \W2{4) 
■) InL 



1 + 51n2 



^vr 



(5.12) 



vr^) 



8vr 



vr^ Li2 I exp 



vr^ 



16vr 



Lis exp 



(5.13) 

At the shifted critical point k[^^ (See (15.71) ). /j^-* takes a scaling form, as it has been pointed 
out by Barber (1974). Notice that at the bulk critical temperature is proportional to 
L^^ In L implying that the finite size-scaling hypothesis (|l.ll) breaks down for the mean 
spherical model with Dirichlet boundary conditions at = 3, while in more realistic models 
the finite-size scaling is valid and the critical Casimir amplitude can be estimated (Krech & 
Dietrich 1992). 

Another expression for the scaling behaviour of the singular part of free energy density 
was obtained in (Barber 1974). In our notations it reads 



1 2 iP' ^ T y"^ 

- Kv m L H 

2 ^ 8vr 2 



K, 



c,3 



\W2{4) 



71n2 



5vr 



+ Qo(y), 



(5.14) 
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where 



Qo{x) 



with 



R{-l)-R 



R{z) 



X 



(5.15) 



In- 



sinh TT 



dw. 



TT\/W 



By comparing (|5.13l) and (15.141) we see that the first three terms in the right hand side 
of (|5.14l) have their counterparts in (15.131) obtained through the replacement y'^ — > — tt^. 
It seems that the terms linear in tt^ and the term proportional to In L were neglected in 
(15.141) . It remains to see what is the situation for the rest of the terms. Unfortunately no direct 
analytic comparison can be made, despite the fact that the integral representation (15.151) can 
be expressed in terms of poly logarithms. To achieve the comparison recourse must be sought 
in numerical methods, so we plot the function Qo{y) and 

1 1 3 1 



D{y) 



IQtt 



(21n2-l) {y'-Tr') 



— A/y2 - Lis fexp 
in \ 



12n 
-2 



1 



IGvr 



Lis exp 



(5.16) 

from (|5.13l) which contains terms that are not present in (15.141) and those that were apparently 
neglected. The result is shown in Figure [TI*! We see clearly that both functions have similar 
behaviours and are shifted one from the other by a constant. This has been checked by 
computing the derivatives of both functions which gives us the same result. The difference 
between the two functions is estimated to be 

AC(3) 



Qo{y)-D{y) 



TT 

16 



-1 + ln(47r2 



(5.17) 



Thus, the scaling functions (|5.13l) and (|5.14l) are equal up to an irrelevant constant that does 
not become singular at the bulk critical point. 

For 3 < < 4, using the asymptotic expansion (|3.21|) . the scaling form of the singular 
part of the free energy density (15.1 II) reads 



+ 



(An] 
2 



(An 



{y'-^'Y 



d 



2-d 



d 



r 



3-d 



°° Kd 

E— 

(iVy 



(5.18) 



Just like the equation for the spherical field we find that the surface contributions to the 
singular part of the free energy violates the standard finite-size scaling hypothesis (11.11) . A 
reformulation of the scaling behaviour would require the introduction of a function with two 

* The numerical evaluation of the named expressions was performed with the aid of WOLFRAM 
MATHEMATICA 6 
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Figure 1. Behaviours of the functions Qo{y) ( 15.14b and D{y) ( 15.16b against y. 



arguments: the scaling variable x and a scaling variable that incorporates the size dependent 
term L'^^^. 



6. Neumann boundary conditions 

6.1. Equation for the spherical field 

To investigate the finite-size effects in the mean spherical model of finite thickness subject 
to Neumann boundary conditions, through equation (12.81 ). we need to estimate the large 
asymptotic behaviour of the sum 

L 



n=l 
^ 2L-1 



+ 2 - 2 COS 



7r(n — 1) 



d-l 



n=0 



■KTl 

+ 2 - 2 COS — 
L - 



-\Wd-i [</' + 4] + iiy,_i 



(6.1) 



The sum in the last line is exactly Sf^^ii'P) (^^^ (13.31) ). Consequently, from (|6.1I) and (13.101) . 
the equation for the spherical field is given by 



K = Wd{<p) - — [Wd-i{<p + 4) - Wd-i 



(47r)'^/2 
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In addition to the restriction on d imposed by the fact that c? = 2 is the lower critical 
dimension and d = Ah the upper one, which is contained in Wd{z) there is another one 
originating from Wd-i{z) that restricts the validity of the asymptotic behaviour (13.121) to the 
interval ?> < d < A. For that reason we will investigate in some details only the cases d = ?> 
and 3 < d < 4. The remaining part of the interval i.e. 2 < < 3 requires special treatment. 

At (i = 3, equation (16.21) takes the simple form 
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'2L 



1^2(4)-— In 2 
47r 
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InL- 
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2Lv/0sinh(Lv/0) 



(6.3) 



A'kL AnL 

This equation was derived in (Dantchev & Brankov 2003) using a method based on (Barber & 
Fisher 1973). At the bulk critical temperature Kc^s the spherical field behaves as ~ L^^ to 
the leading order assuming L^/(J) ^ 1. In the limit L^/(J) ^ 1 one would expect a logarithmic 
behaviour as pointed out by (Dantchev & Brankov 2003). At the shifted critical temperature, 
defined through (Dantchev & Brankov 2003) 
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InL, 



(6.4) 



equation (|6.3I) may be written in a scaling form. An inspection of the expression (16.41) shows 
that the predictions of the theory of finite-size scaling for the shifted critical temperature is 
not fulfilled due to the presence of the term proportional to In L. 

For 3 < d < A the scaling form of the equation for the spherical field (16.21) is given by 



{ATiy/^ 



2 — d\ V^-p f 3 — d 
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(6.5) 



Here again we used the notations x = L~^^^{Kc4 — K) and y = L\f^. As above for the case 
of Dirichlet boundary conditions the standard finite-size scaling hypothesis is violated here 
as well and would need a reformulation in order to take into account the term proportional to 
jj.-'i _ Furthermore, equation (16.51) shows that the leading asymptotic behaviour of the finite- 
size shift of the critical temperature is L^, which is not consistent with the predictions of the 
theory of finite-size scaling. 



6.2. Free energy density 

The free energy density (12.61) in the case of Neumann boundary conditions is obtained by 
analysing the sum 
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(6.6) 
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which contains the term 'P^fl^{(t)), equivalent to (|3.17l) . Using (16.61) and (|3.18l) we obtain 
explicitly 



/3Ff)(T,L; 
where 



m{T-<P) + ^(5F\%{<P) 



1=1 



{21L^) 



(0) = if/,_i M-if/.-i [0 + 4] 



(6.7) 



(6.8) 



is the surfaces contribution. 

Again we are faced with a situation where we have to restrict the validity of 
our expressions due the expansion (13.211) . which should hold for Ud{z) and Ud-i{z) 
simultaneously. Here we consider the cases d = 3 and 3 < d < 4 separately. 

At d = 3, using (|3.21l) and (15.121) the singular part of the free energy density (16.61) reads 
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(6.9) 



Here again we have a logarithmic dependence on L and surface contributions leading to the 
violation of the finite-size scaling hypothesis (11.11 ). This makes the spherical model unsuitable 
for the evaluation of the Casimir amplitude for 0{n) systems with Neumann boundary 
conditions. The value of this quantity is known from renormalization group (Griineberg & 
Diehl 2008). 

For 3 < (i < 4, the singular part of the free energy reads 
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(6.10) 

This suggests that standard finite-size scaling hypothesis (11.11) is not valid for the mean 
spherical model confined to a film geometry with Neumann boundary conditions. Similar 
to the case of Dirichlet boundary conditions a modified finite-size scaling hypothesis is 
necessary to get the appropriate scaling behaviour. 



7. Neumann-Dirichlet boundary conditions 

7.1. Equation for the spherical field 

To extract the finite-size effects from (12.81) in the case of mixed boundary conditions i.e. 
Neumann and Dirichlet on the bounding surfaces, the sum to be considered is 
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2L 
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2 cos 
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, - - (7.1) 

By inspection of the sum in the second line we see that the analysis of this case is tightly 
related to that of antiperiodic boundary conditions. In other words this sum is exactly 



c(") / 



+ 4^^)) with 

2L+1 ) defined in (I4.2al) meaning that in this case the finite- 



size corrections correspond to a film of thickness 2L + 1 subject to antiperiodic boundary 
conditions. Then, along lines similar to the analysis of Section |4l for the equation for the 
spherical field we get, for arbitrary d. 
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(7.2) 



to the leading order in L ^. As before we will restrict ourselves between the lower and the 
upper critical dimensions. For 2 < d < A, using (|3.12l) . we have 
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(7.3) 



With y = L\f^ and x = L^I^{Kc,d — K). The solution of this equation depends explicitly on 
the size L which leads us to the conclusion that the standard finite-size scaling is violated here 
as well. To resolve this issue one would need a modified finite-size scaling assumption. This 
remains valid also for the corresponding finite-size shift of the critical temperature whose 
leading asymptotic behaviour is L^^, rather the than the predicted L^^^^ . Remark however 
that for 0? = 3 the term L'^^^ vanishes and the finite-size scaling is expected to hold. Indeed, 
at d 



47rx 



3, (|7.3I) turns into the simple form 
In 2 cosh 



(7.4) 



This equation is identical to equation (4.31) of (Dantchev & Brankov 2003) obtained by a 
method adapted after (Barber & Fisher 1973). The positive solution of (|7.3I) may be written 
in the scaling form 

y = 9{ND)i>i), (7.5) 
where the universal scaling function is given by 
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At the bulk critical temperature, x = 0, we find the critical amplitude 

yo = 9{ND){^) = 1.45684.... 



(7.7) 



7.2. Free energy density 



The free energy density from equation (12.61) in the case of Neumann-Dirichlet is obtained as 
a result of the analysis of the sum 
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(7.8) 



So apart from a surface term (the last term) the asymptotic form of the free energy has a similar 
expression to (I4.4al) in the case of the antiperiodic boundary conditions with the replacements 
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and L — > 2L i.e. 
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with the surface contribution 
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[Ud{(r) - Ud-i{cr + 4)] 



(7.10) 



For 2 < d < 4, the singular part of the free energy density follows also from that for 
antiperiodic boundary conditions (14.101) and the corrections originating from the surfaces. In 
this case we get 
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(7.11) 

Here again there is a term proportional to L'^^^ that violates the standard finite-size scaling 
hypothesis (|l.ll) in the vicinity of Kc^d and its modification is necessary to get the correct 
scaling. At d = 3 there is no explicit dependence in L. In this case, the critical amplitude 
of the singular part of the free energy density is obtained from (|7.1 II) . Thus the Casimir 
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amplitude for the three dimensional mean spherical film subject to Neumann-Dirichlet 
boundary conditions is 

^(ivo) ^ 0.01922 ... . (7.12) 

This exact result is in conformity with that of (Krech & Dietrich 1992) obtained using 
renormalization group. Similar to the case of antiperiodic boundary conditions the Casimir 
amplitude is positive but smaller in magnitude indicating a weaker repulsing force between 
the surfaces bounding the system. 

8. Discussion 

We investigated the finite-size effects in the d dimensional ferromagnetic mean spherical 
model of finite thickness L subject to different kinds of boundary conditions: periodic (p), 
antiperiodic (a), Dirichlet (D) and Neumann (N) on both surfaces bounding the model, and 
a combination of Neumann and Dirichlet on each surface (ND). We proposed a method for 
the computation of the finite-size corrections of the free energy for arbitrary dimension. Our 
analysis showed that for Dirichlet and Neumann boundary conditions the finite-size effects 
are essentially equivalent to the case of periodic boundary conditions for a film of thickness 
2L and additional surface terms. Similarly, the case of Neumann-Dirichlet was found to be 
related to the case of antiperiodic boundary conditions with thickness 2L. 

The free energy density and the equation for the spherical field were computed for a 
film with arbitrary dimension d subject to the different boundary conditions. In the particular 
case d — 3, our general expressions for (D), (N) and (ND) reduce to those obtained by 
(Barber & Fisher 1973) and (Danchev et al. 1997, Dantchev & Brankov 2003). It is found 
that the singular part of the free energy density has the standard finite-size scaling form for 
2 < d < A only in the cases (p) and (a) i.e. for those boundary conditions which do not break 
the translation invariance of the model. In these cases we estimated the critical amplitude 
of the singular part of the free energy and obtained the values: A^^^^ = — 2^(3)/(57r) = 
—0.153051 . . . and A^") = 0.274543 ... for (p) and (a), respectively. Interpreted in terms 
of the Casimir effect this imply that in the case (p) we have fluctuation-induced attraction 
between the surfaces bounding the model and a repulsion in the case (a). 

For a film subject to (D) or (N) the critical point of the film is shifted from the bulk 
one by surface and finite-size terms, whose leading asymptotic behaviour is proportional to 
L~^. This results disagrees with the scaling behaviour of the finite-size shift predicted by the 
theory of finite-size scaling. In the vicinity of the bulk critical point the standard finite-size 
scaling is not valid in general. At = 3, the solution of the equation for the spherical field 
is proportional to L'^ for (D) and to L^^ for (A^), while the singular part of the free energy 
exhibits a logarithmic dependence on L. For 3 < o? < 4, the standard finite-size scaling for 
the singular part of the free energy needs to be modified. 

In the case (ND), the film has its own critical point shifted from the bulk one by surface 
terms, as well as finite-size effects. The finite-size shift of the critical temperature and the 
scaling form of the free energy do not conform with the theory of finite-size scaling. For 
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2 < d < A, the standard finite-size scaling has to be modified in the neighbourhood of the 
bulk critical temperature. A surprising fact is that the finite-size scaling hypothesis is again 
valid at = 3. At Tc, the square root of the equation for the spherical field is equal to 
1.45684 . . . L"^ and the Casimir amplitude is found to be A^^^) = 0.01922 . . . i.e. a weaker 
than the case (a) fluctuation-induced repulsion. 

It is well known that the mean spherical model is not able to capture the gross features 
of 0{n) models when it is subject to boundary conditions that break the translational 
invariance. To solve this problem remaining in the framework of the spherical model one 
would introduce additional spherical fields to ensure the proper behaviour of the surface spins 
(Singh et al. 1975). Otherwise one can try recovering the equivalence between the spherical 
model and 0{n) spin models by imposing spherical constraints ensuring the same mean 
square value for all spins of the system (Knops 1973). In the case of a film geometry this is 
equivalent to having a spherical constraint on each layer of the system with a space dependent 
spherical field along the finite direction (Hikami & Abe 1976, Bray & Moore 1977, Ohno & 
Okabe 1983) whose relaxed version reduces to the model under consideration. Even in this 
case an accurate study of the problems related to finite-size scaling remains rather untractable 
analytically (see e.g. (Hikami & Abe 1976, Brezin 1983). 
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